Measurement and simulation of the muon-induced neutron yield in lead 



L. Reichhart'' *, A. Lindote'', D.Yu. Akimov^ H.M. Araujo'', E.J. Barnes'' ^ V.A. Belov'', A. Bewick'', A. A. Burenkov^ V. Chepel'', 
A. Cun-ie'', L. DeViveiros'', B. Edwards'^ ^ V. Francis% C. Ghag^*', A. Hollingswol-th^ M. Horn'' ^, G.E. Kalmus^ A.S. Kobyakin^ 

A.G. Kovalenko'', V.A. KudryavtsevS, V.N. Lebedenko^ M.I. Lopes'', R. Liischer'^, R Majewski^ A.St J. Murphy^ F. Neves'", 
S.M. Paling% J. Pinto da Cunha'', R. Preece% J.J. Quenby'', PR. ScovelP^ C. Silva'', V.N. Solovov'', N.J.T. Smith'^'', PE Smith^ 

V.N. Stekhanov^ T.J. Sumner'', C. Thorne'', R.J. Walker''^ 

"School of Physics t-r Astronomy, SUPA University of Edinburgh, UK 
LIP-Coimbra &■ Department of Physics of the University of Coimbra, Portugal 
'^Institute for Theoretical and Experimental Physics, Moscow, Russia 
, ''High Energy Physics Group, Blackett Laboratory, Imperial College London, UK 

I , 'Particle Physics Department, STFC Rutherford Appleton Laboratory, Chilton, UK 

^ ■ ^High Energy Physics Group, Department of Physics &■ Astronomy, University College London, UK 

' ^Department of Physics & Astronomy, University of Sheffield, UK 

OO Abstract 

A measurement is presented of the neutron production rate in lead by high energy cosmic-ray muons at a depth of 2850 m water 
' equivalent (w.e.) and a mean muon energy of 260 GeV The measurement exploits the delayed coincidences between muons and 
^ the radiative capture of induced neutrons in a highly segmented tonne scale plastic scintillator detector Detailed Monte Carlo 
simulations reproduce well the measured capture times and multiplicities and, within the dynamic range of the instrumentation, 
C/2 the spectrum of energy deposits. By comparing measurements with simulations of neutron capture rates a neutron yield in lead 
. ^ of (5.8±0.2) X 10"^ neutrons/muon/(g/cm^) has been obtained. Absolute agreement between simulation and data is of order 25%. 
^ Consequences for deep underground rare event searches are discussed. 
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'1. Introduction 

Rare signal searches, such as those performed for direct dark 
matter detection (^^^ and references therein) and neutrino- 
less double beta decay experiments (f?) and references therein), 
are typically carried out in deep underground laboratories . The 
rock over-burden of such facilities removes or dramatically re- 
duces many of the background signals that would be present 
if the experiments were conducted in surface laboratories. As 
improved sensitivity is achieved, the need to characterise and 
mitigate remaining backgrounds becomes ever more important. 
One of the most problematic backgrounds that still remains is 
that of cosmic-ray muon-induced neutrons, which may become 
a limiting factor in some next-generation rare event searches. 
This specific type of background already shows its impact in 
current dark matter experiments, with XENON 100 reporting 
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it to be the dominant contributor to their nuclear recoil back- 
ground expectation fl. 

Neutrons arising from radioactive decays, for example in 
a fission process or produced in (a,n) reactions following 
a-decay of trace contaminations of heavy radio-isotopes, have 
energies limited to a few MeV. In contrast, neutrons produced 
through interaction of high energy cosmic-ray muons with mat- 
ter can reach energies of several GeV. Consequently, while ra- 
dioactivity neutrons may be effectively controlled by appropri- 
ate shielding constructions and selection of radio-pure build- 
ing materials, removing cosmic -ray induced neutrons is more 
difficult, with the most effective solution being to go deep un- 
derground where the muon flux is reduced by several orders of 
magnitude compared to that at the surface. Eurther mitigation 
of this background involves large muon vetoes, such as instru- 
mented water tanks, to efficiently detect muon tracks far away 
from the detector. 

The neutron production cross-section for high energy muons 
is very large in high- A materials. Yet, several rare event 
search projects utilise large amounts of lead to provide shield- 
ing against ambient y-rays. Thus, the accurate knowledge of the 
production rate of neutrons by cosmic-ray muons in this mate- 
rial is very important for assessing and planning the capability 
of these projects, present and future. 

A penetrating cosmic-ray muon may produce neutrons via 
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four main processes: (i) muon spallation — muon-nuclear inter- 
action via the exchange of a virtual photon, resulting in nuclear 
disintegration, (ii) muon capture (only dominant for shallow 
depths, <100 m w.e.), (iii) photo-nuclear interactions in muon- 
triggered electromagnetic showers, and (iv) hadron-production 
in hadronic cascades initiated by the muon. These secondary 
cascades make up most of the muon-induced neutron produc- 
tion in deep sites. Specifically, neutrons are predominantly cre- 
ated by photo-nuclear interactions of y-rays produced in elec- 
tromagnetic showers, neutron inelastic scattering, pion spalla- 
tion and pion absorption at rest. The rate of neutron production 
by direct muon nuclear interaction is significantly smaller than 
for the other processes listed [6,- lOJ. 

The non-trivial task of measuring the cosmic -ray muon in- 
duced neutron yield has been pursued by a number of under- 
ground experiments (see lH - IT | for a compilation of such re- 
sults). Most recently, the KamLAND collaboration has pre- 
sented muon-induced neutron rates for a number of target 
isotopes I12il . Additional work from other groups is ongo- 
ing i^. While for low- A targets agreement between the dif- 
ferent measurements and simulation toolkits (GEANT4 lll3ll . 
FLUKA 114 ,15.1 ) is reasonable, studies of heavy targets are 
somewhat controversial and inconsistent iITgIi . Older mea- 
surements for Pb targets |17 lH, including beam measure- 
ments [19], without Monte Carlo simulations of neutron pro- 
duction, transport and detection, show much larger neutron 
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yields than expected from simulations |8- 10, 20]. On the con- 
trary, measurements with the veto of the ZEPLIN-II experi- 
ment at the Boulby Underground Laboratory showed an over- 
production in the simulation by ~80% (if]. 

Here we present a new measurement of the muon-induced 
neutron yield in lead using the data accrued by a highly 
segmented anti-coincidence detector installed around the 
ZEPLIN-III dark matter instrument. The measurement was 
conducted in parallel to the 319-day long second science run 
of the experiment in 2010/1 1. 

2. Experimental apparatus 



The ZEPLIN-III instrument 11221 12311 is a dual phase liquid/gas 
xenon detector, built to observe low energy nuclear recoils re- 
sulting from elastic scattering of weakly interacting massive 
particles (WIMPs). The final scientific exploitation of the long- 
running ZEPLIN project at the Boulby Underground Labora- 
tory (at 2850 m w.e.) achieved cross-section limits for scalar 
and spin-dependent WIMP-neutron channels of 3.9 x 10"^ pb 
and 8.0 x 10"^ pb near 50 GeV/c^ (90% confidence), respec- 
tively 11241 12511 . For the second science run the detector was 



upgraded with a new array of photomultiplier tubes (PMTs), 
decreasing internal background sources significantly 112 6ll , and 
a 1 -tonne plastic scintillator anti-coincidence detector system 
(the 'veto'), mounted around the main instrument [27]. De- 
signed for rejecting background events in the WIMP target, the 
veto also helped to decrease systematic uncertainties in the es- 
timation of background rates due to independent measurements 
of y-ray and neutron rates in the vicinity of the dark matter 
detector ll28ll . In addition, the veto provided an independent 
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Figure 1 : (Colour online) CAD rendering of the veto system surrounding the 
ZEPLIN-III dark matter detection instrument. The veto barrel consists of 32 
vertical Gd-loaded polypropylene pieces (white) suiTounded by the same num- 
ber of active scintillator modules (black), with PMTs housed in cups and re- 
cessed into the lower polypropylene structure. The roof of the veto detector is 
composed of 20 scintillator modules, which are placed on top of a roof plug. 
The lower polypropylene structure contains no Gd and rests on a copper and 
lead base. Finally, a lead castle (only the first few lead blocks are shown on 
the sides facing the back) envelops the entire assembly (~2.3 m in length and 
~2.4 m in height). For display purposes only, a quarter of the scintillator bars 
from the barrel are not drawn to reveal the ZEPLIN-III detector 



measurement of the high energy cosmic-ray muon flux and its 
spallation products. The performance of the instrument is well 
understood through data and validated Monte Carlo simula- 
tions izt 



The veto system consists of 52 individual plastic scintilla- 
tor modules (51 were active during the second science run) 
surrounding a 15 cm thick Gd-loaded polypropylene shield- 
ing, which encircles the ZEPLIN-III instrument. This entire 
structure is then enclosed in a 20 cm thick lead castle. A CAD 
rendering of the full setup is shown in Fig. [1] For a detailed 
description of the design and performance of each individual 
component see Ref. ll27ll . Here, only a short summary is pre- 
sented. 

The structure formed by assembling the individual modules 
can be described by two main geometrical shapes: a circu- 
lar barrel composed of 32 vertical scintillator bars and a roof 
constructed from 20 individual scintillator blocks. Each barrel 
bar has a trapezoidal cross-section with parallel sides of length 
15 cm and 12 cm and a height of 15 cm. The length of the 
barrel scintillators is 1 m. The roof sections are of four differ- 
ent lengths (80, 75, 67 and 51 cm) oriented to form a pseudo- 
circular shape divided into quadrants and are of rectangular 
cross-section with side lengths of 15 cm x 16 cm. The indi- 
vidual detector bars are made from polystyrene-based plastic 
scintiUator UPS-923A (p-terphenyl 2%, POPOP 0.02%), pro- 
duced by Amcrys-H, Kharkov, Ukraine ll29il . A single PMT 
(ETEL-9302 KB) is optically coupled to one end of each indi- 
vidual scintillator bar. Additionally, all bars have been wrapped 
in PTFE sheet of high diffuse reflectivity, and a highly-specular 
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reflective aluminised Mylar film is located at the end opposite 
to the PMT to increase light collection. 

The polypropylene shielding inside the scintillator construc- 
tion is loaded with ~Q.4% Gd by weight |28]. Thus, many neu- 
trons, moderated to thermal energies, undergo radiative cap- 
ture on '''^Gd due to its very high capture cross-section of 
2.4 X 10^ barn HQ]. This is of great advantage for detecting and 
identifying radioactivity neutrons from internal detector com- 
ponents. These are detected with high efliciency through the 
emission of 3^4_y-rays (totalling ~8 MeV) with a mean delay 
of only ~1 1 jUS 1281] . For the more energetic muon-induced neu- 
trons, which are mostly produced externally, a slower capture 
on hydrogen in the plastic scintillator is expected. 

Data were accrued with a dedicated data acquisition system 
(CAEN model V1724), digitising waveforms with 14-bit res- 
olution, an input range of +2.25 V, 40 MHz bandwidth and a 
sampling rate of 10 MS/s. The waveforms of recorded events 
were 320 fj.s in length; they were parameterised using a bespoke 
data reduction software ('RaVen') adapted from that developed 
for the ZEPLIN-III instrument B^ . 

The veto detector was operated in 'slave' and 'master' mode 
simultaneously. In slave mode the veto acquisition system was 
triggered by an external signal generated by ZEPLIN-III. The 
trigger point and timeline lengths were tailored to enable quasi 
dead time free recording of coincident events. The master mode 
allowed for independent triggering of the veto system when cer- 
tain requirements were met. One of these conditions was the 
sum of simultaneously occurring pulses in the roof modules ex- 
ceeding a set threshold (summed in a dedicated hardware unit). 
At this depth, most cosmic-ray muons have an arrival direction 
which is close to vertical, and thus such a trigger condition pro- 
vides a high efficiency for detection of cosmic-ray muons, but 
adds little to the total data storage or rate implied for the exper- 
iment. 

Critical for a long running experiment is the stability of the 
detector system over time. Thus, a number of parameters, in- 
cluding electronic gains (measured with the single photoelec- 
tron response of the PMTs), coincidence rates, background 
rates, tagging efficiencies of electron recoil events and environ- 
mental parameters, were monitored throughout the course of 
the experiment. Additionally, a dedicated calibration run was 
performed on a weekly basis, with a pulsed blue LED, coupled 
via fibre optic cable to each individual scintillator bar at the end 
opposite to the PMT. Monitoring of the mean of the single pho- 
toelectron peak and of the centroid of the LED-generated peak 
over the duration of the experiment confirmed the system's sta- 
bihty Jli. 



3. Monte Carlo simulations 

Simulated primary muon energy spectra and angular distribu- 
tions were obtained by propagation of atmospheric muons from 
the Earth's surface through an appropriate depth of rock using 
the MUSIC code [32, 331; this distribution was then sampled 
with the MUSUN code [sstl • The energy, momentum, posi- 
tion and charge of each muon was recorded at the point where it 



intersected the surface of a cuboid fully enclosing the main cav- 
ern of the laboratory. The cuboid included an extra 5 m of rock 
on each side, except for the top which enclosed a total of 7 m 
of additional rock. The mean energy of the muon distribution 
was ~260 GeV and 20 million of these muons were generated. 
The equivalent live-time of the final simulation for the present 
study amounts to ~3.1 years. 

The comprehensive simulation that was developed for the 
ZEPLIN-III experiment has already been well established in 
previous studies |23, 2^ 34 1. Complementary investigations 



of the veto detector have also been performed This 
simulation was updated to run with version 9.5 (patch 01) of 
GEANT4 for this work. 

To model the physical processes for this setup the modular 
physics list Shielding, currently recommended for shielding 
applications at high energies, was implemented. It uses the 
Fritiof string model (FTP) and the Bertini cascade (BERT) for 
the high and low energy ranges (up to 5 GeV), respectively, 
similar to the FTFP_BERT reference list but with different neu- 



tron cross-section data (JENDL-HE-2007 |135|] up to 3 GeV 
and evaluated cross-sections [36] above 3 GeV) iIItIi . Neu- 
tron interactions below 20 MeV are described by high-precision 
data-driven models with data obtained from the ENDF/B-VII 
library (3^. Additionally, thermal scattering off chemically 
bound atoms was implemented for neutron energies below 4 e V, 
which is especially important to model thermalisation in the 
plastics 

Secondary particle production thresholds ('cuts') were set 
to 0.1 mm for y-rays and e"/e^ which, in lead, translate to 
~30 keV and ~250 keV, respectively. This is safely below 
photo- and electro-nuclear reaction thresholds. 

The output generated by the simulation has been designed 
to recreate that of the experiment, i.e. a waveform-like read- 
out with a resolution of 0.1 //s for all 52 individual channels 
separately. Thus, direct comparison to data as well as the use 
of similar analysis cuts for experimental and simulated data is 
possible. 

4. Event selection 

During the second science run, it was required that the veto be 
maximally sensitive to the low energy deposits expected from 
multiply scattering radioactivity neutrons and y-rays. Conse- 
quently, bias voltages for each PMT were adjusted to deliver 
a dynamic range in the region of 1-70 photoelectrons (phe), 
corresponding to approximately 20-1300 keV at the far end of 
the scintillator A minimum-ionising muon crossing the full 
thickness of a scintillator bar deposits at least ~20 MeV and, 
thus, muon signals, and a greater number of MeV energy de- 
posits from ambient y-ray background, result in heavily satu- 
rated pulses. Given a single range data acquisition, recording 
of non-saturated muon events simultaneously with the signal 
expected from captured neutrons would not be possible. Se- 
lection of muons from this data set is therefore non-trivial, but 
can be achieved by searching for coincident saturated signals in 
roof and barrel scintillators, due to the optical separation of the 
modules. 
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Highest energy deposition in a roof module [phe] 



Figure 2: (Colour online) Highest energy deposition observed in one of the roof 
scintillator modules versus the time-coincident highest energy deposition in a 
barrel scintillator module. The dashed line indicates the graphical cut used to 
select the muons. Note that energy depositions >100 phe are saturated. 



Figure |2] shows the greatest energy deposition observed in 
a roof module plotted against the largest corresponding (co- 
incident) signal in a barrel module for each event, occurring 
within ±0.2 fis around the trigger point. This is similar to the 
prompt coincidence window used for tagging y-ray events in 
ZEPLIN-III 1 28]. A well separated population is observed, 
with a graphical selection criterion indicated. Here, the pho- 
toelectron scale is defined by using a constant conversion factor 
between the pulse area and the pulse height parameter, which 
was utilised for the single photoelectron calibration. Given that 
the pulse area is less affected by saturation than the pulse height 
(due to the abrupt cutoff in the latter), the impact of saturation 
can be pushed to higher energies (>100 phe), and so improve 
separation of event populations. 

Selection of muons in the simulation followed a very similar 
procedure. Firstly, events with a minimum energy deposition 
observed from the summed signal of the veto roof, analogous 
to the trigger function of the veto detector, were selected. Ad- 
ditionally, as in the data, a cut on time coincidence (0-0.4 jus) 
between roof and barrel was applied. In Fig. |3] the Monte Carlo 
data are plotted as a function of the largest energy deposition in 
the (coincident) barrel module only. Separate curves are shown 
for all events satisfying the coincidence condition, and for only 
those events corresponding to energy depositions directly re- 
sulting from muon traversal of scintillator modules. The differ- 
ence between the two curves is predominantly due to the energy 
depositions from particles generated in showers as muons pass 
nearby. A simple cut at the position indicated by the dashed 
line selects a population which is composed of ~93% muon 
energy depositions with the required coincidence, i.e at least 
one roof and one barrel module firing within the defined coin- 
cidence window. 

Confirmation that the identified region in the experimental 
data corresponds to the muon event region in the Monte Carlo 
is provided by comparing the event distributions between pairs 
of roof modules and barrel modules, as shown in Fig. |4] Here, 
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Highest energy deposition in coincident barrel module [phe] 

Figure 3: (Colour online) The plot shows the highest energy deposition ob- 
served in a barrel scintillator module when measured in coincidence with a 
roof module for all prompt energy depositions in the simulation (black solid). 
The red dashed spectrum shows the same but for muon hits only, ('. e. the muon 
crosses both the roof and the coincident barrel module. The cut used to select 
muon events is indicated by the thick dashed vertical line. ~889'o of muons with 
a roof - baiTel coincidence have energies above this threshold. 



the two upper panels show the distribution of roof modules 
(numbered 32-5 1 ) registering a coincidence with a specific bar- 
rel module (modules 3 and 19, as indicated). Similarly, the 
lower panels show which barrel modules are in coincidence 
with which roof modules (39 and 46). The combination of the 
relative orientations of the modules with respect to each other, 
their individual response functions, and the asymmetric impact 
of the surrounding laboratory geometry, results in a complex 
distribution of coincidences between modules. However, the 
Monte Carlo reproduces the experimental data reasonably well 
(the average reduced value of these 5 1 coincidence contri- 
butions is ~ 1 .9), confirming that the selected experimental data 
correspond to cosmic -ray muon events. 

An overall efficiency for pure muon events, including pre- 
viously mentioned effects and the (geometric) requirement for 
coincidence between barrel and roof is 36.8+0.6%, where the 
error includes uncertainties due to the precise choice of the lo- 
cation of the selection cuts. 

A total number of 7979 muons was selected from the 
full dataset translating to a rate of 32.3+0.4 muons/day. 
By comparing the measured rate to the Monte Carlo pre- 
diction, using the normalised flux through a sphere in the 
simulation in a similar way to ll2lll . we deduce a muon 
flux of (3.75+0.09) X 10"*^ muons/s/cml This result is 
in excellent agreement with the last reported value for 
the muon flux in the Boulby Underground Laboratory of 
(3.79+0.15) X 10"** muons/s/cm^ [21], measured in the cav- 
ern hosting both the ZEPLIN-II and ZEPLIN-III detectors, and 
~8% lower than the value obtained for another cavern in Boulby 
reported in ll40ll . 
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Figure 4: (Colour online) Sample of coincident channels from two barrel slabs 
(top, module 3 and 19) and two roof slabs (bottom, module 39 and 46) with 
all modules (channels) of the roof and the barrel, respectively. The simulation, 
scaled to the total muon rate observed in the data, is shown by the red dashed 
hatched histogram in comparison to the data (black histogram). 



5. Muon-induced neutron yield 

The vast majority of detected neutrons produced by muons in 
this set-up originates in the ~60-tonne lead shield, which pro- 
tects the experiment from ambient y-rays. To determine the 
muon-induced neutron yield in lead from the present data we 
count the number of neutrons captured in the veto following a 
recorded muon event. This is compared with simulations per- 
formed using the same analysis cuts. We note that a data set 
with single photoelectron resolution is a real asset: at the ex- 
pense of a small increase in background rate, the low threshold 
analysis increases the number of detected neutrons substantially 
in comparison with previous works. 

5.7. Experiment 

As described previously, neutrons are identified through signals 
occurring in one or more of the 5 1 scintillators as a result of the 
y-rays emitted following their capture. These signals are de- 
layed relative to the muon's passage due to the time for thermal- 
isation and capture to occur Ideally, the data would be searched 
for the signatures of neutron captures over the entire period in 
which these signals may arrive. However, the PMT response 
to a large energy deposition is such that the timelines become 
at first heavily saturated, and then are followed by a large sig- 
nal overshoot. For extreme energy depositions the overshoots 
persisted for up to 40 jUS. The effect of these 'dead' waveform 
periods can be seen in Fig. |5] showing a significantly reduced 
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Figure 5: (Colour online) Time delay distribution for all recorded pulses above 
given thresholds relative to the muon. 



pulse rate for the first -40 /is after the muon trigger. Thus, the 
timeline for detecting delayed neutrons was restricted to the re- 
gion of 40-300 yus relative to the observed muon. An efficiency 
of ~47% was retained from this timeline selection cut (calcu- 
lated from simulations). Furthermore, the maximum number 
of recorded pulses was restricted to 300 entries per event (an 
equivalent cut was implemented in the analysis of the simula- 
tion). The impact of this restriction is discussed in Section lS!2l 

For the detection of muon-induced neutrons, as compared to 
internally-generated radioactivity neutrons, which are also of 
concern in dark matter searches, one expects an increased im- 
portance of neutrons capturing on hydrogen. In the dark matter 
veto case, the neutrons have scattered within the ZEPLIN-III 
device and thus have a high geometrical likelihood of being 
captured in the Gd-loaded polypropylene shielding immedi- 
ately surrounding the target. Most muon-induced neutrons 
come from outside of the setup and will more likely be cap- 
tured in the hydrocarbon scintillator material surrounding the 
Gd-loaded shielding. A single ~2.2 MeV y-ray is emitted fol- 
lowing capture on hydrogen, and therefore signals observed in a 
single scintillator module are more likely to occur - in contrast 
to the several y-ray signature from Gd capture, which can be 
recorded simultaneously in several scintillator modules. Due 
to the relatively long capture times of neutrons on hydrogen 
in comparison to captures on gadolinium, the rejection of the 
first 40 fis of the waveforms reduces the probability of detec- 
tion with single scintillator signals by only -26%. 

Single scintillator events are more exposed to backgrounds, 
and careful consideration of thresholds and a good knowledge 
of those backgrounds are required. Since available statistics of 
the limited pre-trigger timeline fraction are very scarce, an addi- 
tional data set from the same run with similar trigger conditions 
was used to estimate the background correctly. A dataset of 
synchronised (with ZEPLIN-III) 'slave' triggered veto events 
(see Section |2]i was considered to calculate the background. 
Optimisation of the number of neutron captures observed in 
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Figure 6: (Colour online) Time delay distribution of channel multiplicity one 
events (threshold of ^10 phe) in ZEPLlN-111 coincident background data. The 
dashed (red) line indicates the constant fit to the background. 



Table 1 : Measured number of neutrons per muon from the data in comparison 
to neutron rates extracted from simulations using the same requirements and 
cuts as in the experimental data analysis. Background rates, for correction of 
the data, are listed individually for the different channel multiplicities, with their 
required thi'esholds detailed in the text. The errors given for the data are the sum 
of statistical eiTors and the rate coming from random accidental coincidences 
of pulses calculated from the average observed pulse rate for a given threshold. 
EiTors of simulated rates are statistical only. 







Data 




Simulation 


Channel 


Events/ 


Background n/muon 


n/muon 


mult. 


muon 


rate 


(bkg.corr) 




1 


0.216(5) 


0.019(1) 


0.197(5) 


0.145(2) 


2 


0.088(3) 


0.0049(5) 


0.083(3) 


0.076(1) 


3 


0.039(2) 


0.0019(3) 


0.037(2) 


0.0321(9) 


>4 


0.029(2) 


0.0008(2) 


0.028(2) 


0.0231(8) 


Total 


0.372(7) 


0.026(1) 


0.346(7) 


0.275(3) 



the muon triggered data, with respect to the number of false 
events due to background, results in a threshold of ^10 phe 
being chosen for single scintillator events (the present results 
were shown to be largely insensitive to the precise threshold). 
Figure |6] shows the rate of background events with a threshold 
of ^10 phe applied. The rate is approximately constant, i.e. it 
is independent of the time since the trigger occurred. 

Following the methodology used in the analysis of the dark 
matter search data 12511 , coincident signals in multiple scintil- 
lators can also be searched for, detecting multiple scatters and 
y-rays following neutron capture on gadolinium at later times. 
Coincidences are defined as occurring within +0.2 fis of each 
other To optimise efficiency, different signal size thresholds 
have been required depending on the number of signals in coin- 
cidence, balanced against the rate of false signals arising from 
non-neutron related sources (background and induced noise). 
Noise from the PMTs can be intrinsic, i.e. from thermionic 
emission and internal radioactive decays, or directly induced. 
Especially after larger signals, such as resulting from muons, 
positive ions generated from ionisation of residual gases in the 
PMTs lead to secondary signals, creating afterpulses at short 
time scales of up to several fis dependent on the ion transit time 
(see ll4lll42ll and references therein). In the present data after- 
pulses of small amplitudes are suppressed at short decay times 
due to the large pulse overshoots observed following a muon en- 
ergy deposit (see Fig.|5j. A second noise component observed 
at longer time scales, visible in Fig. |5]between ~30-50 fis af- 
ter the start of the muon signal, with sizes of 4 phe and be- 
low, may be attributed to the organic scintillator Luminescence 
with long time constants is expected from phosphorescence and 
dela yed fluorescence processes in the plastic scintillator (see 



e.g. iillillii). 

These additional signals could lead to false coincidences 
between scintillators, generating spurious neutron detections. 
Based on the event rates, the probability of false coincidences 
can be calculated. It was found that for neutron capture events 
with a channel multiplicity of two, i.e. two scintillator bars fir- 



ing within +0.2 jUS of each other, a signal size requirement of 
threshold ^4 phe (in each pulse) was sufficient to remove af- 
terpulses. For three-fold coincidences between scintillators, a 
threshold of ^2 phe per signal was found to be appropriate, 
and for four or more scintillators, a threshold at the level of a 
single photoelectron was sufficient. For consistency, a global 
requirement was set that regardless of the number of scintilla- 
tors fired in coincidence, all events must have a total signal size 
of at least 8 phe. Despite the lower threshold for multiple scin- 
tillator events, accidental rates arising from background are a 
lot smaller due to the required coincidence of pulses. The same 
dataset used earlier to estimate the background rate in the single 
scintillator case has thence been utilised to calculate the contri- 
bution from background to the yields of neutron captures found 
from the multiple scintillator requirements. Background rates 
found are at the level of statistical uncertainties. 

Table [1] summarises the results. Each instance in which the 
designated criteria were met is interpreted as indicating a neu- 
tron capture. Most muon-induced neutron captures are ob- 
served through events seen in single scintillators only, despite 
the higher threshold required. However, a significant num- 
ber also generate energy depositions observed in coincidence 
in several scintillators. Overall, a mean of 0.346+0.007 neu- 
trons (including background corrections) are observed for every 
muon detected. 

5.2. Comparison with simulations 

Monte Carlo simulations of the experiment have been per- 
formed as described in Section |3] The dimensions and param- 
eters of th e appar atus have been previously measured and doc- 
umented 11271 |34j, |45tl, including signal gains and attenuation 
lengths of the scintillators so that photoelectron spectra can be 
generated. This allows Monte Carlo pseudo-data to be analysed 
using identical routines as used for the real data. The overall 
agreement on the rate of detected neutrons between data and 
simulation obtained in this work is good, at the level of 25%. 
For the initial discussion the total number of neutrons has been 
normalised to the data. 
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Figure 7: (Colour online) Time delay distribution of detected captured neutrons 
from experimental data (black solid histogram) and simulations (red dashed 
hatched histogram). The constant background has been subtracted from the 
data histogram. Results from simulation are normalised to the total number of 
neutrons observed in the data. 



Figure 9: (Colour online) Energy depositions of detected neutrons from back- 
ground corrected data (black solid) and simulations (red dashed hatched his- 
togram) below the saturation point in the data, i.e. the energy scale is given in 
absolute number of photoelectrons (1 phe =: 20 keV). Results from simulation 
are normalised to the total number of neutrons observed in the data. 
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Figure 8: (Colour online) Comparison of channel (scintillator module) multi- 
plicities per detected neutron in the data (black solid) to simulations (red dashed 
hatched histogram). The data are background corrected according to Table [T] 
Results from simulation are normalised to the total number of neutrons ob- 
served in the data. 



Figure 10: (Colour online) Relative fraction of neutron multiplicities per muon, 
i.e. the number of delayed signals observed after a muon trigger in the defined 
time window, for background corrected data (black sohd) and simulation (red 
dashed hatched histogram) nonnalised to the total number of observed muons 
in each case. 



Figure |7] compares the time delay distributions for detected 
neutrons from data (solid black) and simulation (red dashed). 
Excellent agreement between the two distributions is found. 
Moreover, the module (channel) multiplicity per neutron event 
can be used as an additional consistency check, beyond the ini- 
tial muon identification, taking advantage of the segmented na- 
ture of the detector Figure|8]shows the number of channels with 
coincident signals involved in each individual neutron event. 
The data are corrected for the contributions from background 
coincidences, as given in Table [1] Again, excellent agreement, 
over the full range of channel multiplicities, is demonstrated. 

In Fig.|9]the energy depositions associated with the observed 
(captured) neutrons are given in the region before the onset of 
saturation, with excellent agreement between simulation and 
data obtained. Although the energy calibration is only known 
to within 10% due to, amongst other factors, the saturation 



of the data (see also previous studies with the same instru- 
ment 1 28, 45J), tests in varying the energy scale by this amount 
resulted in only small neutron rate differences and are consid- 
ered in the systematic error of the simulated rate. 

A muon may produce more than one fast neutron in a cas- 
cade, resulting in several neutron capture signals at different 
times and in different locations in the veto. Figure [10] shows 
the relative fraction of observed neutrons per muon for data and 
simulation. When exploring neutron multiplicities, rather than 
scaling the simulation to the total number of neutrons observed 
in the data, a simple normalisation to the number of detected 
muons has been applied. Background corrections assume even 
distribution of background events. As such, most non-neutron 
signals occur in one of the empty waveforms following a muon 
trigger, making up almost 90% of all observed muon events. 
Generally, good agreement is observed. 
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When scaled to the number of neutrons detected, the simu- 
lations reproduce well the time distributions, the energy depo- 
sitions and the number of scintillators involved in each event. 
The absolute numbers of neutrons expected to be observed per 
muon, as determined from the simulation for each individual 
channel multiplicity, are also summarised in Table [T] show- 
ing an overall reduced neutron rate from simulations of ~20% 
{i.e. the total yield from the data exhibits an increase over the 
simulation by ~26%). Discrepancies are largest for single scin- 
tillator events. At higher multiplicities absolute agreement be- 
tween simulation and data is of order 10-20% (cf. ~36% for 
single scintillator events). 

The expected total muon-induced neutron rate calculated 
from simulations is 0.275+0.003 (stat.) +°JJo7 (syst.) neu- 
trons/muon. Systematic errors are calculated from the variabil- 
ity in the energy calibration. 

As previously mentioned, a restriction on the maximum 
number of recorded pulses was applied to the data, and simi- 
larly in the analysis of the simulation. Importantly, the number 
of selected muon events, in both data and simulation, affected 
by this limitation are in excellent agreement with each other, 
further supporting the performance of the Monte Carlo simula- 
tion. When including all energy depositions in the simulation a 
higher absolute neutron rate is observed (~30%). This increase 
can be associated with only a few muon events (approximately 
1 in 1700) featuring exceptionally high neutron multiplicities. 
It is worth noting that these high multiplicity events are less sig- 
nificant for dark matter searches due to the generally high veto 
efficiencies expected for these. 

Table |2] shows the relative production of neutrons in differ- 
ent materials for all neutrons generated in the simulation and 
for detected neutrons only. As expected, nearly all neutrons 
are produced in the rock cavern of the underground laboratory, 
reflecting that the simulation included sufficient volume to re- 
move edge effects. Importantly, less than 1.5% of detected neu- 
trons are produced in the rock, confirming the effectiveness of 
the shielding setup of the ZEPLIN-III detector On the other 
hand, the lead component of the shielding enclosure provides 
an effective target for neutron production by high energy comic- 
ray muons, with ~95% of all neutrons created there. 

Table [3] lists the specific elements involved in the capture of 
the neutrons, both for detected neutrons and for all the neutrons 
in the simulation. As previously mentioned, the vast majority of 
detected muon-induced neutrons are captured on hydrogen, em- 
phasising the importance of measuring the single -2.2 MeV y- 
ray from this process. The captures on Gd amount to ~7.0% for 
this configuration at these neutron energies (being much more 
effective for internal radioactivity neutrons due to the detector 
geometry described in Section|2}. 

5.3. Muon-induced neutron yield in lead 
As shown in Table |2] the detected neutrons have predominantly 
been produced in lead. Thus, the observed neutron rate may 
be used to derive an absolute neutron production yield in this 
material. The methodology used follows that of Refs. ¥k Elll . 
and is essentially to scale an idealised simulation of neutron 
production by a mono-energetic beam of muons in pure lead 



Table 2: Fractions of neutrons produced in different materials for all generated 
neutrons in the simulation and for detected neutrons only. 



Production material of 

Material all neutrons detected neutrons 

Lead 02% 951)% 

Rock 99.8% 1.4% 

Steel - 1.2% 

CgHg - 0.9% 

Copper - 0.8% 

CH2 - 0.5% 

Gd-epoxy - 0.1% 

Liquid Xe - 0.1% 



Table 3: Fractions of neutrons captured on different elements for all and for 
detected neutrons only. 

Capture element of 

Element all neutrons detected neutrons 

'- 711% 

Fe - 11.5% 

CI 94% 7.0% 

Gd - 7.0% 

Pb - 1.3% 

C - 1.1% 

Cu - 0.6% 

Na 6% 0.2% 

Mn - 0.2% 



by the ratio in rate observed between the present data and the 
full detector simulation (assuming that the fraction of detected 
neutrons produced in lead (~95%) is well described by the sim- 
ulation). 

The simulation of a mono-energetic muon beam in lead 
was conducted as follows. Neutron production was recorded 
for a mono-energetic 260 GeV ji' beam (mean muon en- 
ergy at Boulby), incident on the centre of a lead block of 
3200 g/cm^ thickness. Figure [TT] shows the differential en- 
ergy spectrum of neutrons produced. Only neutrons from 
the central half length of the lead block were considered to 
avoid surface/edge effects. To prevent double counting of 
neutrons in neutron inelastic processes, the first neutron pro- 
duced in each reaction was dismissed independent of its en- 
ergy. As summarised in Table |4] a neutron production rate 
of (4.594+0.004) x lO"-' neutrons/muon/(g/cm^) was obtained 
for the physics list and version of GEANT4 used through- 
out this work (Shielding with version 9.5). However, since 
the experimental muon-induced neutron rate was found to 
be a factor of 1.26+0.03 (stat.) +°°3 (syst.) higher, our re- 
sults suggest a true production rate by 260 GeV muons of 
(5.8+0.1 (stat.) j (syst.)) x 10"-' neutrons/muon/(g/cm^), as- 
suming neutron transport and detection are modelled accu- 
rately. 

In this analysis, uncorrected arrival of muons is assumed, 
as opposed to muon bundles produced together by primary 
cosmic -rays in the atmosphere. A study based on a simple ap- 



8 




10-^ 10' 1 10 10^ 10' 10^ 
Neutron energy [MeV] 



Figure 11: Differential energy spectrum of muon-induced neutrons produced in 
lead from fi" of 260 GeV. 



proximation to find the survival probability of muons at a given 
depth showed the effect to be negligible and the error to be very 
small for the measurement performed with the ZEPLIN-III 
veto detector 

A similar comparison between simulation and experiment 
was performed for the ZEPLIN-II anti-coincidence system JtI 
I21I1 . In that work a muon-induced neutron yield in lead of 
(1.31+0.06) X 10"^ neutrons/muon/(g/cm^) was reported. Here 
we have revisited the simulation, now using GEANT4 version 
9.5 and the Shielding physics list, including thermal scatter- 
ing cross-sections, and a significantly larger sample of primary 
muons. This resulted in a new estimate for the neutron yield in 
lead of (3.4+0.1) x 10"^ neutrons/muon/(g/cm^) in that setup. 
While it is clear that a significant contribution of the newly ob- 
tained ZEPLIN-II yield comes from the updated simulation, 
there remains a significant discrepancy with the present result. 
One possible explanation for this is that the angular distribu- 
tion of emitted neutrons may not be accurately modelled. The 
virtually complete geometrical enclosure of the ZEPLIN-III 
veto with respect to the lead shielding would be expected to 
make the present results less susceptible to such effects than in 
ZEPLIN-II. 

We have also explored the evolution of the neutron produc- 
tion yield with successive versions of GEANT4. To do this fur- 
ther simulations of a mono-energetic /i -beam focused on a lead 
block have been performed. Table |4]summarises the results, in- 
cluding the yield obtained with version 8.2 from Ref. ||7|. In 
addition, combination of different physics lists and GEANT4 
versions are listed, also linking the custom list used in [7] to the 
current high energy reference lists. The bespoke physics list is 
very similar to QGSP_BIC_HP, featuring the Quark-Gluon String 
(QGS) theoretical model at high energies coupled to nuclear de- 
excitation with a pre-compound model, the intra-nuclear Binary 
Cascade (BIC) model below 6 GeV and a data-driven high- 
precision neutron package (NeutronHP) to transport neutrons 
below 20 MeV down to thermal energies. Reasonable varia- 
tion of change-over energies between the BIC and QGS models 



Table 4: Muon-induced production yields for neutrons for different versions 
of GEANT4 and physics lists (for 260 GeV muons). The neutron yield from 
version 8.2 is based on the value reported in Ref. |7]. A small modification 
has been applied to correct for a previously unaccounted error in the rejection 
of neutrons produced in neutron inelastic processes to avoid double counting 
(refen'ed to as 'stars' in that work). 

GEANT4 physics list muon-induced neutron yield 



version [neutrons/muon/(g/cm^)] 

8^2 custom built (2.846+0.006) x 10"^ 

9.4 custom buih (3.304+0.003) x 10"^ 

9.4 QGSP_BIC_HP (3.376+0.003) X 10"^ 

9.4 Shielding (3.682+0.003) x 10"^ 

9.5 QGSP_BIC_HP (3.993+0.004) x 10"^ 
9.5 QGSP_BERT_HP (4.369+0.004) x 10"^ 
9.5 FTFP_BERT (4.467+0.004) X 10"^ 
9.5 Shielding (4.594+0.004) x 10'^ 



in the custom physics list in comparison to the reference one 
have little impact (<3%) on the overall neutron yield. A steady 
increase with every new version of GEANT4 is demonstrated. 

The Shielding physics list shows not only the largest 
muon-induced neutron production yield in comparison to other 
reference lists, but is also subject to the highest increase in go- 
ing from version 9.4 to 9.5. This is explored in detail in Fig. [12] 
showing the individual contributions from the most important 
neutron creation processes for muons in lead. The main differ- 
ence lies in the increased neutron production in inelastic scat- 
tering of hadrons and in particular neutrons. A ~38% higher 
production yield for this process is observed. 

Part of the increase observed between versions 9.4 and 9.5 
of the toolkit (applicable to all standard lists used in this 
study) can be attributed to the muon-nucleus interaction model 
(G4VDMuonNuclearModel); as in previous versions, this still 
relies on the Kokoulin mu-nuclear cross-sections i46ll . but the 
final state of the hadronic vertex is now replaced by a 7r° 
interacting further through the Bertini intra-nuclear cascade. 
The previous model (G4MuNuclearInteraction) replaced 
the virtual photon with n'^^^ instead, which would then interact 
through the low/high energy parameterised models (LEP/HEP) 
— these are known to yield fewer neutrons. There has also been 
increased neutron production in the FTF model, which may ac- 
count for some of the enhanced yields in the Shielding and 
FTFP_BERT lists; The addition of the Reggeon cascade 
which can cause more nucleon secondaries, is a possible ex- 
planation, but further study is required [,48.1 . 

6. Conclusion 

For the development of future rare-event searches, especially 
in the context of direct dark matter experiments, accurate data 
on muon-induced neutron yields in several materials is of great 
importance, as is the ability to simulate these processes using 
modern Monte Carlo toolkits. Complex models inform the de- 
sign of large and expensive shielding and veto systems around 
these experiments, as well as the interpretation of their data 
(background expectations). There exists significant uncertainty 
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Figure 12: (Colour online) Absolute neutron yields of the most important pro- 
duction processes for muon-induced neutrons generated from filing 260 GeV 
fi^ on lead using the Shielding physics list and GEANT4 version 9.5 (black 
histogram) and version 9.4 (red dashed histogram). The neutron creation pro- 
cesses are: 1: photo-nuclear interaction of y-rays (y —> N), 2: neutron inelastic 
scattering (n — » N), 3: pion spallation (n —> N), 4: muon spallation (ji — »N), 5: 
proton spallation (p N), 6: pion absorption (n^ abs) and 7: all other neutron 
production processes. 



in the simulated muon-induced neutron rate, as evidenced by 
the steady change in the total neutron yield with every new ver- 
sion of GEANT4 and physics list and available data have been 
inconclusive. 

In this study, a dataset from 319 days of operation of the 
ZEPLlN-111 anti-coincidence detector has been analysed for 
high energy cosmic -ray muons. The number of muon-induced 
neutrons has been evaluated by detecting delayed y-ray signals 
following radiative captures. A muon flux in the Boulby Un- 
derground Laboratory of (3.75 ±0.09) x 10 muons/s/cm^ has 
been determined, consistent with and improving upon previ- 
ous measurements. The muon-induced neutron detection rate 
was measured to be 0.346+0.007 neutrons/muon traversing the 
ZEPLlN-111 scintillator veto. Monte Carlo simulations, using 
GEANT4 (version 9.5) and the Shielding physics list with 
the same cuts and thresholds applied as used for the analysis of 
the data, resulted in a neutron capture rate of 0.275;^^ ^^^ neu- 
trons/muon (quadratically combined statistical and systematic 
errors), which is ~2Q% lower than the experimentally measured 
value. However, absolute rates aside, the simulation reproduced 
very well all tested parameters, strengthening confidence in the 
results. The ratio of neutron rates between data and simula- 
tion have been used to evaluate a muon-induced neutron yield 
in pure lead of (5.8+0.2) x lO"-' neutrons/muon/(g/cm^) for a 
mean muon energy of 260 GeV. Additional simulations explor- 
ing previous versions of the GEANT4 simulation package con- 
firm the trend of an increasing neutron production rate in lead 
with every successive distribution of GEANT4 (also shown in 
other simulation studies 49]). A further increase by 25% is 
encouraged by the results presented in this work. 

Finally, as our results confirm, the very significant contribu- 
tion of lead in the production of muon-induced neutrons, the use 



of lead based shielding to prevent y-rays from the environment 
to propagate into the sensitive volume of the detector should be 
carefully assessed for any future rare event search. Alternative 
shielding compositions, such as large water tanks surrounding 
the detectors, are already being used in some current dark mat- 
ter searches [50-52], as well as discussed for near future next 
generation experiments ]j53 , 54 ] . 



Acknowledgement. The UK groups acknowledge the 
support of the Science & Technology Facilities Coun- 
cil (STFC) for the ZEPLlN-111 project (in particular, the 
analysis work presented here was supported by STFC grants 
ST/K006436/1, ST/K003178/1, ST/K003208/1, ST/K006428/1, 
ST/K003 186/1, ST/K006444/1 and ST/K006770/1) and for 
maintenance and operation of the underground Palmer Lab- 
oratory which is hosted by Cleveland Potash Ltd (CPL) at 
Boulby Mine, UK. The project would not have been possible 
without the co-operation of the management and staff of CPL. 
Additionally, we want to thank the Boulby science facility team 
for their support during underground aspects of this work. We 
also acknowledge support from a Joint International Project 
award, held at ITEP and Imperial College, from the Russian 
Foundation of Basic Research (08-02-91851 KO.a) and the 
Royal Society. LlP-Coimbra acknowledges financial support 
from Funda9ao para a Ciencia e Tecnologia (FCT) through 
the project-grants CERN/FP/109320/2009 and /116374/2010, 
and postdoctoral grants SFRH/BPD/27054/2006, /47320/2008, 
/63096/2009 and /73676/2010. Furthermore, we acknowledge 
the Edinburgh Compute and Data Facilities for accommodating 
the heavy usage of the Edinburgh computer cluster 'Eddie' for 
this analysis. This work was supported in part by SC Rosatom, 
contract contract #H.4e.45.90.11.1059 from 10.03.2011. The 
University of Edinburgh is a charitable body, registered in 
Scotland, with registration number SC005336. 



References 

[1] D. Akimov, Techniques and results for the direct detection of dark matter 
(review). Nuclear Instruments and Methods in Physics Research A 628 
(2011) 50-58. 

[2] E. Figueroa-Feliciano, Direct detection searches for WIMP dark matter. 
Progress in Particle and Nuclear Physics 66 (201 1) 661-673. 

[3] V. Chepel, H. Araujo, Liquid noble gas detectors for low energy particle 
physics, arXiv: 1207.2292 [physics. ins-det]. 

[4] A. Barabash. Double Beta Decay Experiments, Physics of particles and 
nuclei 42 (2011) 613-627. 

[5] E. Aprile, et al., Dark Matter Result from 225 Live Days of XENONIOO 
Data, Physical Review Letters 109 (2012) 181301. 

[6] V. Kozlov, et al., A detection system to measure muon-induced neutrons 
for direct dark matter searches, Astroparticle Physics 34 (2010) 97-105. 

[7] A. Lindote, et al., Simulation of neutrons produced by high-energy muons 
underground, Astroparticle Physics 31 (2009) 366-375. 

[8] V. Kudryavtsev, N. Spooner, J. McMillan, Simulations of muon-induced 
neutron flux at large depths underground, Nuclear Instruments and Meth- 
ods in Physics Research A 505 (2003) 688-698. 

[9] H. Araiijo, V. Kudryavtsev, N. Spooner, T. Sumner, Muon-induced neu- 
tron production and detection with GEANT4 and FLUKA, Nuclear In- 
struments and Methods in Physics Research A 545 (2005) 398^1 1. 



10 



[10] Y.-F. Wang, et al., Predicting neutron production from cosmic-ray muons, 
Pliysical Review D 64 (2001) 013012. 

[11] D.-Mei, A. Hime, Muon-induced background study for underground lab- 
oratories, Physcial Review D 73 (2006) 053004. 

[12] S. Abe, et al.. Production of radioactive isotopes through cosmic muon 
spallation in KamLAND, Physical Review C 81 (2010) 025807. 

[13] S. Agostinelli, et al., GEANT4- a simulation toolkit. Nuclear Instruments 
and Methods in Physics Research A 506 (2003) 250-303. 

[14] G. Battistoni, et al., The FLUKA code: Description and benchmarking, 
Proceedings of the Hadronic Shower Simulation Workshop 2006, Fer- 
milab 6-8 September 2006, M. Albrow, R. Raja eds., AIP Conference 
Proceeding 896 (2007) 31^9. 

[15] A. Ferrari, P. Sala, A. Fasso, J. Ranft, FLUKA: a multi-particle transport 
code, CERN-2005-10 (2005), INFN/TC.05/11, SLAC-R-773. 

[16] V. Kudryavtsev, L. Pandola, V. Tomasello, Neutron- and muon-induced 
background in underground physics experiments. The European Physical 
Journal A 36 (2008) 171-180. 

[17] G. Gorshkov, V. Zyabkin, R. Yakovlev, Production of neutrons in Pb, 
Cd, Fe and Al by high-energy muons, Soviet Journal of Nuclear Physics- 
USSR 18 (1974) 57-61. 

[18] L. Bergamasco, S. Costa, P. Picchi, Experimental results on neutron pro- 
duction by muons at 4300 m we, Nuovo Cimento della Societa Italiana di 
Fisica A-Nuclei Particles and Fields A 13 (1973) 403-412. 

[19] V. Chazal, et al., Investigations of fast neutron production by 190 GeV/c 
muon interactions on different targets. Nuclear Instruments and Methods 
in Physics Research A 490 (2002) 334-343. 

[20] M. Marino, et al.. Validation of spallation neutron production and prop- 
agation within GEANT4, Nuclear Instruments and Methods in Physics 
Research A 582 (2007) 611-620. 

[21] H. Araujo, et al.. Measurements of neutrons produced by high-energy 
muons at the Boulby Underground Laboratory, Astroparticle Physics 29 
(2008)471^81. 

[22] D. Y. Akimov, et al.. The ZEPLIN-m dark matter detector: Instrument 
design, manufacture and comissioning, Astroparticle Physics 27 (2007) 
46-60. 

[23] H. Araujo, et al.. The ZEPLIN-IH dark matter detector: Performance 
study using an end-to-end simulation tool, Astroparticle Physics 26 

(2006) 140-153. 

[24] V. N. Lebedenko, et al.. Results from the first science run of the ZEPLIN- 
111 dark matter search experiment. Physical Review D 80 (5) (2009) 
052010. 

[25] D. Y. Akimov, et al., WIMP-nucleon cross-section results from the second 
science run of ZEPLlN-111, Physics Letters B 709 (2012) 14-20. 

[26] H. M. Araujo, et al., Radioactivity backgrounds in ZEPLIN-in, Astropar- 
ticle Physics 35 (2012) 495-502. 

[27] D. Y. Akimov, et al.. The ZEPLIN-Ill anti-coincidence veto detector, As- 
troparticle Physics 34 (2010) 151-163. 

[28] C. Ghag, et al.. Performance of the veto detector incorporated into the 
ZEPLIN-m experiment, Astroparticle Physics 35 (2011) 76-86. 

[29] A.-H. Corp, www.amcrys-h.com. 

[30] M. Chadwick, et al., ENDF/B-VII.O: Next generation evaluated nuclear 
data library for nuclear science and technology, Nuclear Data Sheets 
107 (12) (2006) 2931-3059. 

[31] F. Neves, et al., ZEBRA: the ZEPLIN-III Reduction and Analysis pack- 
age. Journal of Instrumentation 6 (2011) PI 1004. 

[32] P. Antonioli, et al., A three-dimensional code for muon propagation 
through the rock: MUSIC, Astroparticle Physics 7 (1997) 357-368. 

[33] V. Kudryavtsev, Muon simulation codes MUSIC and MUSUN for un- 
derground physics. Computer Physics Communications 180 (2009) 339- 
346. 

[34] M. Horn, et al, Nuclear recoil scintillation and ionisation yield in liquid 
xenon from ZEPLIN-IH data, Physics Letters B 705 (201 1) 471-476. 

[35] JENDL High Energy File 2007, http://wwwndc.jaea.go.jp/ffipnd/jendl/jendl- 
he-2007.html. 

[36] V. Barashenkov, Reprint: Nucleon - Nucleus Cross-sections, Joint Insti- 
tute for Nuclear Research (1990) P2-89-770. 

[37] T. Koi, A GEANT4 Physics List for Shielding Calculations, at the 
SATIFIO at CERN (June 2010). 

[38] Cross Section Evaluation Working Group, ENDF/B-Vn data files may be 
accessed through the national nuclear data center www.nndc.bnl.gov. 

[39] S. Gamy, et al., GEANT4 Transport Calculations for Neutrons and Pho- 



tons Below 15 MeV, IEEE Transactions on nuclear science 56 (4) (2009) 
2392-2396. 

[40] M. Robinson, et al., Measurments of muon flux at 1070 m vertical depth 
in the Boulby underground laboratory. Nuclear Instruments and Methods 
in Physics Research A 511 (2003) 347-353. 

[41] J. Birks, The Theory and Practice of ScintUlation Counting, Pergamon 
Press Limited, 1964. 

[42] L. Campbell, Afterpulse measurement and correction. Review of Scien- 
tific Instimients 63 (1992) 5794-5798. 

[43] U. Kawada, J. Ito, Q.-W. Wang, Temperature dependence of spurious 
pulses in use of plastic scintillation detectors. Applied Radiation and Iso- 
topes 60 (2004) 403^07. 

[44] T. Marvin, Single-Photoelectron Noise Reduction in Scintillation De- 
tectors, IEEE NUCLEAR SCIENCE SYMPOSIUM AND MEDICAL 
IMAGING CONFERENCE RECORD 1-3 (1996) 764-767. 

[45] L. Reichhart, et al., Quenching factor for low energy nuclear recoils in a 
plastic scintiUator, Physical Review C 85 (2012) 065801. 

[46] A. Bogdanov, et al., Geant4 Simulation of Production and Interaction of 
Muons, IEEE Transactions on nuclear science 53 (2) (2006) 513-519. 

[47] V. Uzhinsky, Development of the Fritiof Model in Geant4, Joint Interna- 
tional Conference on Supercomputing in Nuclear Applications and Monte 
Carlo 2010 (SNA + MC201()), Tokyo, Japan. 

[48] D. Wright (SLAG, GEANT4 Collaboration), private communication. 

[49] V. Tomasello, Background Simulations for a Large-scale Cryogenic Dark 
Matter Experiment, Ph.D. thesis, University of Sheffield, High Energy 
Physics Group (November 2009). 

[50] D. Akerib, et al.. The Large Underground Xenon (LUX) Experiment, Nu- 
clear Instruments and Methods in Physics Research A 704 (2013) 111- 
126. 

[51] M. G Boulay for the DEAP Collaboration, DEAP-3600 Dark Matter 
Search at SNOLAB, Journal of Physics: Conference Series 375 (2012) 
012027. 

[52] H. Sekiya, for the XMASS collaboration, XMASS, Journal of Physics: 

Conference Series 308 (2011) 012011. 
[53] D. MalUng, et al.. After LUX: The LZ Program, arXiv:1110.0103v2 

[astro-ph.IM]. 

[54] H. Kraus, EURECA - the European Future of Dark Matter Searches with 
Cryogenic Detectors, Nuclear Physics B (Proc. Suppl.) 173 (2007) 168- 
171. 



11 



